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Abstract 

We study various constraints in the non-decoupling MSSM scenario. The light Higgs boson may evade 
the direct search experiments at LEPII or Tevatron while the 125 GeV Higgs-like boson is identified as 
the heavy Higgs boson in the spectrum. Two direct consequences of the scenario are the flavor violation 
induced by the light charged scalar and the spin-independent scattering between neutralino and nuclei in 
dark matter direct detection experiments. With combined flavor constraints B — > X s "f and B s — > \i + \T~ and 
direct constraints on Higgs properties, we find best fit scenarios with light stop of 0(500 GeV), negative A t 
around -750 GeV and large /i-term of 2-3 TeV However, large parameter region in the survival space under 
all bounds may be further constrained by H — » rr if no excess of rr is confirmed at LHC. We only identify 
a small parameter region with significant H — >■ hh decay that is consistent with all bounds and reduced rr 
decay branching fraction. In addition, if current dark matter mostly consists of neutralino, direct detection 
experiments like XENON 100 also puts stringent bound over this scenario with light Higgs bosons. The 
light stops which are required by flavor constraints can further enhance the scattering cross section. 
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I. INTRODUCTION 



A Higgs-like boson of 125 GeV has been discovered at the Large Hadron Collider (LHC) at 
CERN via two cleanest channels, the di-photon (gg — > h — > 77) and four-lepton (gg — > h — > 



ZZ* 



with i, j = e ± , / u. ± ) modesUJ]. Later both ATLAS and CMS collaboration also 



■ % ~ % ^ 



reported observations in the di-lepton (gg — > h — > WW* — > ifvitjvj) channels with the mass 
range consistent with the four-lepton measurement [2L However, the confirmation of whether it 
is the Higgs boson of the standard model (SM) will require comprehensive and precise measure- 
ments of Higgs properties. The deviation of the Higgs couplings from the SM ones may imply 
the existence of the beyond SM physics, in particular, the excess in the di-photon channel with 
0"obs./VsM ~ 1.5 — 2.0 at both ATLAS and CMS. In all extension theories, additional charged and 
neutral scalars are inevitable. Therefore, searches of other Higgs-like states also provide direct 
test to models beyond SM physics. The LHC and Tevatron collaborations Q3D have put stringent 
bounds over the SM Higgs, particularly heavy Higgs decaying into pure leptonic final states via 
WW and ZZ. For instance, CMS collaboration has excluded the SM Higgs of 110-121.5 GeV 
and 128-600 GeV at 95% C.L. The LEPII experiments also exclude the SM Higgs with mass lower 
than 114 GeV via e + e~ — > Zh channel. These bounds at the same time apply to various models 
with Higgs extension. 



For two decades, weak scale supersymmetry has been the most elegant candidate to cancel the 
quadratic divergence if the Higgs boson is indeed a fundamental scalar. Within the supersymmetric 
framework, there exist several scenarios where the di-photon decay branching fraction is enhanced, 
for instance, models with light stau|4J] or light stop qO. Another particularly interesting region of 
non-decoupling limit in minimal supersymmetric standard model (MSSM) has been discussed by 
various authors Q6Wl4£| . It was observed that there might exist even lighter Higgs h which evades 
the search at LEP [ Jdue to suppressed ZZh coupling and thus production of Zh. The light Higgs 
h can then have Mh < mz while the Higgs-like boson of 125 GeV can be identified as the heavier 
degree of freedom H . To reduce the ZZh coupling gzzh = sin(/3 — a) which is the vacuum 
expectation value (vev) of h, simple realization is to let h be the H d -\ike boson since large m t 
naturally requires large v u . Given h is a mixture state as — sina(Re Hd) + cosa(Re H u ), this 
scenario prefers sin a ~ — 1 and large tan (3 which suppresses the v d . In the limit of large tan /? as 
sin f3 — > 1, sin a — > — 1 gives the sin(/3 — a) approaches zero. On the other hand, within MSSM, 



at tree level, the Higgs mass matrix gives 

tan 2a M\ + m% 



(1) 



tan 2/3 M\ - m| 

Taking Ma — >■ and the (3 — > n/2 as limit of large tan (3, one can get a — > —n/2 which results in 
sin(/3 — a) — > and reproduce the previous requirement. For M A > 200 GeV, the gzzh goes to the 
SM value. However, since the charged Higgs state H ± are at the similar scale as M# as M% ± = 
M\ + at tree level, small Ma leads to lighter H ± which suffers from direct search bounds 
of light charged Higgs. Drell-Yan production of charged Hig gs p air at LEP e + e _ — > H + H~ put 
strict bounds as M H ± > 80 ~ 100 GeV depends on its decay fl 1 5M . Combining all the constraints, 
one expect an intermediate Ma region around mz scale to be consistent with the LEPII Zbb search 
and charged Higgs search at LEP, Tevatron and LHC. In the limit of Ma — > m z , h and H masses 
at tree level are degenerate which is known as non-decoupling limit. 

By requiring Mh to be at 125 GeV, the first consequence of these non-decoupling scenarios 
is that the charged Higgs is around similar scale. Charged scalar below top quark mass receives 
stringent bound from the ATLAS search of t — > bH + with H + — >• t + u requires the BR(t — > 
bH + ) x BR(#+ t + v t ) < 1 ~ 5% for mass range M H ± in between 90 and 160 Gevfl]. 
In the conventional Two-Higgs-Doublet models (2HDM) such a light charged Higgs suffer severe 



constraints due to flavor violation processes lll7ll . For example, one might be concerned by B u — >• 
tv t and B — > D^ru T decays which receive charged Higgs contributions at the tree-level. The 
two most sensitive parameters involved in Higgs interaction are Ma and tan (3. As we argued, Ma 
is taken to be not much heavier than m z and LEP2 Zh search prefers a relatively large tan (3. In 
addition, as we will show later, the recent search of t — > bH + at the LHC restricts tan/3 ~ 10 
in non-decoupling region. For B u — > rv T decay, The W^-mediated SM contribution is helicity 
suppressed. Therefore, even though the charged scalar is somewhat heavier, its contribution could 



be comparable to the SM part if tan /3 is not small J18|-|20j]: 



BR(£?+ -> t+z/)mssm - 2 x 2 



BR(£?+ -> r+v) SM V M% + ^) (2) 
where the MSSM corrections to the down quark and lepton mass matrix have been neglected, 
which is safe for tan (3 ~ 10. For M H + lies around 120 ~ 150 GeV, the MSSM prediction would 
be about 20% ~ 30% smaller than the SM result of (0.95 ± 0.27) x 10~ 4 . While the experimental 



world average is (1.65 ± 0.34) x 10 4 before 2012 112 ill . Belle updated their measurement at 



ICHEP2012 with much smaller value 0.72+^? x 10~ 4 for hadronic tag of r UM- So in the non 



decoupling limit, a light charged Higgs with tan (3 ~ 10 is well consistent with the new Belle 



measurement. Similarly, the charged Higgs contribution to B — > D^ru T decays are not very 
significant in the interesting region of M H + and tan (3. Therefore we will not discuss the bounds 
from B + — >■ t + u and B — > D^*>tv t decays further in our study. 

On the other hand, the penguin b — > s processes are also sensitive to the charged Higgs effects. 
Generally, b — > 57 and B s — > (i + fJ,~ are two most stringent constraints. But choosing appropriate 
MSSM parameters, supersymmetric contributions may cancel part of the SM and charged Higgs 
amplitudes |23|]. For example, b — >■ s transition mediated by the scalar top quark (stop) loop in 
MSSM may cancel the top quark loop in SM and 2HDM in some parameter region. One can thus 
expect that light stop in MSSM may significantly reduce the flavor violation 112411 . In this paper, 
we start with this argument and study whether scenarios with light stop can resolve the tension in 
flavor physics due to the light charged Higgs H ± . 

Search of /i — > &y at the MEG experiment will soon reach BR( / u — > ej) ~ 1 x 1CT 13 . The 
one loop contribution from charged state to — > is suppressed by small lepton masses and 
additional helicity-flip. The largest contribution in Higgs mediated ji — > is usually the Barr- 
Zee two-loop effects involving the charged scalar coupling to a top-bottom loop. However, I125n 
shown that the charged Higgs contribution only reach the sensitivity for tan of 60 for M A of 
100 GeV where the tan (3 is much larger than what is considered in non-decoupling scenarios. 

With conserved R-parity, the thermal relic abundance of the lightest neutralino (LSP) can often 
be identified with dark matter (DM), consistent with the current cosmological observations. In 
recent years, direct detection of weakly interacting (WIMP) DM particle through the DM scat- 
tering with nuclei has excluded large parameter space of supersymmetric DM and put stringent 
bound on many models. The latest bound from XENON100 is about 5 x 1CT 9 pb for DM mass 
around 200 GeV Q26Q. Neutral Higgs states h, H can also mediate the scattering between DM 
and nuclei which is of 1/M^ H . Then the second consequence of non-decoupling scenarios is that 
the spin-independent scattering is significantly enhanced by the interaction through neutral Higgs 
H, A of O(100 GeV) MM- Therefore, models with only neutralino DM in the non-decoupling 
MSSM suffer stringent constraints from direct detection experiments. In addition, light stop which 
may significantly improve the flavor physics behavior of non-decoupling MSSM as argued above, 
would further enhance the scattering of DM and nuclei and put stronger bound on non-decoupling 
scenarios with only neutralino DM^. 

In the next section, we discuss some general constraints on the non-decoupling scenarios and 
1 If the DM is not dominated by the neutralino component, the bound can be evaded. 
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the scan results. Then we discuss in details the physics interpretation of the scan results, in par- 
ticular, light stop contribution to cancel light charged Higgs and its implication to Mh, di-photon, 
di-tau decay and the direct detection experiments of neutralino dark matter. We then conclude in 
the final section. 

H. GENERAL CONSTRAINTS AND THEIR IMPLICATIONS 

In this section, we first scan the parameter space with focus on non-decoupling region with Ma 
is at 0(mz) then discuss in details the physics interpretation of scan results. 

Latest data from the LHC require the resonance to be at 125 GeV with di-photon decay en- 
hanced with respect to the SM prediction. We therefore impose the selection rules as 

• M H : 125 ± 2 GeV; 

• ^77 = ^Zbs/^SM : 1 ~ 2; 

• Combined direct search bounds from HiggsBound3.8.0; 

• BR(fi -»■ X sl ) < 5.5 x 10" 4 ; 

• BR(£ S -> n+ ir) < 6 x 1(T 9 . 

Without loss of generality, we fix masses of the following sfermions as 

M h = M& a = Ms = Mr = M e - , = 1 TeV , (3) 

Ql,2 U1 >2 "1,2,3 -^1,2,3 e l,2,3 ' V > 

and the gauginos as 

Mi = 200 GeV, M 2 = 400 GeV, M 3 = 1200 GeV . (4) 

As argued, our study focus on the flavor constraints of the non-decoupling MSSM and b — > s 
transitions like B — > X s ^ and B s — > [i+ provide the most severe constraints. Light stop usually 
helps to cancel the charged Higgs contribution in b — > s transition. On the other hand, for light 
stop below 500 GeV, we find that the gluon fusion production of H is suppressed significantly with 
respect to the SM value due to the cancellation between top squark and top quark in the loop. Thus, 
for light stop {Mi < 500 GeV), it is difficult to achieve enhanced di-photon. For comparison, we 
take the third generation up quark masses as 

Mq 3 = M- t = 500 GeV and a second group with 1 TeV . (5) 
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We do the scan over four parameters I 

Ma ■ 95 ~ 150 GeV 
tan/3 : 1 ~ 30 

(j, : 200 GeV ~ 3 TeV 
A u = A d = A e : -3 ~ 3 TeV . (6) 

Discussed by many authors[|4|], light stau states may significantly enhance the di-photon rate of 
the Higgs-like boson decay which are observed by both ATLAS and CMS collaborations. On 
the other hand, we don't require much stronger di-photon bound as i? 77 to be the experimental 
preferred central value of 1.5. Stau states are irrelevant to the flavor constraints from b — > s 
transition but only give minor change to the Higgs boson mass. Therefore, we don't take light stau 
in the study. 

We use FeynHiggs 2.9.2 [Q] 3 with HiggsBounds 3.8.0 H and SUSY_Flavor 2.01 [30] to 
perform the scan here. Figured] shows the scan results in 2D-plot of A t and /i. Ma and tan (3 are 
also varied but aren't shown in the figures. Figure [D-(a) is the heavy stop scenario with Mq 3 = 
Mi = 1 TeV and (b) is the light stop scenario with Mq 3 = M 4 ~ = 500 GeV. Points in red region 
pass the direct search bounds from HiggsBounds with a heavy CP-even Higgs M H = 125 ± 2 
GeV and an enhanced diphoton rate 1 < i? 77 < 2. The points in blue region pass in addition the 
constraint of BR(£> — > X s y), while the points in black region pass all the constraints, including 
further the restriction of BR(£> S — > 

The scenario with heavy stop can survive the B — > X s ^ constraints. However, none of the 
scanned points can pass the B s — > ji + ji^ . In the case of light stop of 500 GeV, we find a small 
survival parameter region with negative A t around 750 GeV and large yU-term between 2 to 3 TeV. 
In the following subsections, we discuss in details the physics implications of the scanned results. 



A. b — > and B s — > fi + fi 

b — > S7 and B s — > turns out to be the most stringent flavor physics bounds in the non- 

decoupling limit. The helicity for the involved quark states must be flipped in b — > 57. Hence, 

2 We confine ourselves to Ma < 150 GeV for larger splitting between h and H which can reduce the tt decay 

branching ratio. Details is discussed later. 

3 In this scan, we take the pole mass of m t instead of the running m t mass. The survival parameter region after scan 

may be shifted by a few percent. 
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FIG. 1. Scan Results in [A t , /z] plane. The heavy (light) stop scenario with Mq^ = = 1 (0.5) TeV is 
shown in the left (right) plot. The red region pass the direct search bounds from HiggsBounds with a heavy 
CP-even Higgs Mh = 125 ± 2 GeV and an enhanced diphoton rate 1 < i? 77 < 2. The blue region pass in 
addition the constraint of BR(f? — > X s ^), while the black region pass all the constraints, including further 
the restriction of BR(B S — > fi + fi^). 

both chiral symmetry U (3)q x U (3)d and electroweak symmetry SU(2) L x U(l)y must be broken. 
The SM b — > s transition is mediated by the charged weak boson W~ and only left-handed quarks 
are involved in the weak interaction. Consequently, b — > 37 is suppressed by mass insertion of 
bottom quark mass nib m the SM. In MSSM, the charged Higgs i/~-top quark loop contribution to 
b — > S7 is also suppressed by m b insertion, and has the same sign as the SM amplitude. Besides the 
above contributions, squarks can also generate b — > S7 in MSSM which may not flip the helicity 
of the involved quark states, for instance, loops with right-handed stop-Higgsino (t^ — H u ) or 
left-handed stop-Wino (t L — W). Therefore the squark contributions, in particular the top squark 
ones, are not necessarily suppressed by m b , which is helpful to cancel the SM and charged Higgs 
amplitudes with appropriate MSSM parameters. Consequently, scalar top quark with small Mj, 
say ~ 500 GeV, could significantly reduce b — > s transition. 

The squark contributions can be decomposed into chargino penguins, wino penguins and gluino 
penguins. Chargino penguins contain tan /3-enhanced term which arises from v u insertion in 
Qd c (H*). The term explicitly breaks Peccei-Quinn symmetry as well as i?-symmetry and is pro- 
portional to [xA t . This contribution would destructively interfere with the SM and charged Higgs 
amplitudes in case of uA t < OlL l32ll. In our study, gluino penguins are also important as they 
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contain terms enhanced by /xtan (3 and terms chirally enhanced by nig /nib. Numerically, we use 
the FeynHiggs program to get the Non-MFV result of BR(£> — > X s j). The experimental world 



average of this process is (3.43 ± 0.22) x 10" 4 Q21 



2_while the SM prediction up to NNLO per- 



turbative QCD corrections is (3.15 ± 0.23) x 10~ 4 fl33h . However, B — > X s ^ decay is evaluated 
only at NLO in the FeynHiggs program, which produces the SM result as 3.8 x 10~ 4 . This is 
about 30% larger than the NNLO SM prediction. Taking this and the theoretical and experimental 
uncertainties into account, we require loosely BR(£> — > X s ^)mssm < 5.5 x 10~ 4 as the selection 
rule in the scan. 



In the SM, BR(£> S — > is strongly helicity suppressed by the small muon mass as 

m Ai/ m B 3 ' wn i cn leads to a tiny branching ratio of (3.27 ± 0.23) x 10~ 9 [|34f | . However, it is 
well known that the MSSM contributions to this decay could be enhanced several orders of mag- 
nitude larger than the SM prediction in large tan (3 limit, as the leading contribution of Higgs 
penguin diagrams to the branching ratio are proportional to tan 6 f5. In our study, tan /3 ~ 10 is not 
very large, so all the 1-loop diagrams have to be considered, including the charged Higgs diagrams 
which is enhanced up to tan 2 at the amplitude level. Notice that B s — > decay is even more 

sensitive to the MSSM parameters in the non-decoupling limit as the neutral Higgs bosons are all 
light. Experimentally, a combined search of ATLAS, CMS and LHCb has set the upper limit of 



370 . this upper limit 



4.2 x 10" 9 Q35J] for time integrated branching ratio. As pointed out in ||36L 
should be reduced by about 10% when compared with the theoretical calculation. Numerically, 
we use the SUSY_FLAVOR program [30J] to get the complete NLO result of BR(B S /i + /i~). 
However, we notice that SUSY_FLAVOR evaluates this branching ratio to be 4.8 x 10~ 9 in the SM. 
This is about 50% larger than the SM prediction of (3.27 ± 0.23) x 10" 9 in Q24D, probably mainly 
due to different choice of hadronic parameters. Taking this into account, we set the corresponding 
selection rule to be BR(i? s — > /i + ^~)mssm < 6 x 10" 9 in the scan. 



In Fig. [Q the black region which satisfy all the constraints give 10 4 BR(£? — > Xs^mssm in 
the region [4.9, 5.3] and 10 9 BR( J B S ii + h~)mssm in the region [2.3, 4.3]. Notice that BR(fi 
X s -f) is always larger than the SM prediction, which is mainly due to the enhancement of light 
charged Higgs. For BR(£? S — > / u + /i~), it is always somewhat smaller than the SM prediction. 
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B. Higgs mass and its decay properties 



We discuss the mass spectrum of the Higgs bosons in non-decoupling MSSM and its decay 



properties in this section. More general discussion can be found in 113 811 . In particular, we focus on 
the parameter region that minimizes the flavor violation in b — > s transition. Combined constraints 
from B — > Xgj and B s — > we take light stop of Mf ~ 500 GeV with negative A t of 

O(-750 GeV) and large //-term of 2-3 TeV. MSSM contains two SU{2) L doublets H u and H d with 
the ratio of their vevs tan (3 = v u j v d . To evade LEPII bounds, non-decoupling limit corresponds to 
a region of much lighter Hd state with small vev in the spectrum. After spontaneously electroweak 
symmetry breaking, MSSM gives rise to five physical states of Higgs bosons, the two CP even 
scalar h, H with one CP odd scalar state A and charged scalars H ± . The two CP even scalar 
bosons h, H arise from mixing of the real gauge eigenstates (Re Ha, Re H u ), 



— sm a cos a 
cos a sin a 




(7) 



After diagonalizing the general mass matrix of neutral Higgs 



M = 2 , (8) 



the masses of two CP-even Higgs are 



Ml = M\ x sin 2 a + A^^cos 2 a - M\ 2 sin 2a, 
Mjj = Ml 1 cos 2 a + Ml 2 sm 2 a + Mf 2 sm2a, 



(9) 



To illustrate the feature, we take the limit of sin(/3 — a) — > which is the vanishing limit of gzzh 
to completely suppress the Zh production at LEPII. As a result of sin a — > — 1 and sin /3 — > 1, we 
have 

-Mil 

(10) 

- M22 



Radiative corrections to the elements in mass matrix Eq. [9] are given in Il39ll . We list the most 
relevant M.22 in Eq. [TT] 

3 



t 



h |^!l2sin 2 /3. 



+ 



167T 2 

A t a 
M 2 

IVI SUSY 



1 + 

A t a 



16tt 2 



12Mg USY 



(1.5^ + 0.5^-8^) 
t 



1 + 



16vr 2 



(3y 



16^) 



v2 Vb ■ 2 o H 

o sin p 

16tt 2 ^ 



M 4 

1V1 SUSY 



t 



1 + ^ (% 2 - %t - 16^ 2 ) 



0(y 2 t m 



(11) 



where #3 is the QCD running coupling constant, y t and yb are the top and bottom Yukawa cou- 
plings. Msusy is the arithmetic mean of top squark masses M^. A t is the SUSY breaking A-term 
associated with top squark and yt, is the Higgsino mass parameter, t is defined as In {Mg USY /m 2 ) 
and 

a = A t — fi/ tan/3 . (12) 

In EqdH only the leading terms in powers of yb and tan (3 have been retained. Even though the 
EqUUis only valid in the limit of small splittings between the running stop masses, it shows the 
qualitative feature for how couplings to stop and sbottom modify the Higgs masses. M4 is cos (3 
dependent which is suppressed in the limit of large tan (3. Therefore, the M H is not very sensitive 
to Ma and with tan (3 ~ 10, varying Ma by 100 GeV results in 10 GeV difference in Mh- Unlike 
the m™ ax scenario with a = \/6M S usy which is usually used in many studies, minimization of 
the flavor violating b — > s transition leads to our best fit parameter region around 



.4, 



-750 GeV, Mi ~ 500 GeV, n ~ 2000 - 3000 GeV . 



(13) 



The particular choices of A t and M t - significantly modifies the Higgs boson masses through ra- 
diative corrections. In our studies, we use the FeynHiggs program to compute the mass spectrum 
of Higgs in which full radiative corrections of Higgs masses have been implemented [28y. Figure 
[2] show how the M h H ^ H ± vary with respect to Ma for one of our benchmark points tan (3 = 11, 
M- t = 500 GeV, A t = -740 GeV and /i = 2300 GeV. For a large range of M A , M H is around 
125 GeV. Non-decoupling limit of nearly degenerate h, H lies near Ma ~ 160 GeV. 

Since H is mostly H u with large tan (3, the v u dominates the electroweak symmetry breaking v. 
The couplings between H and W + W~ and top quark t are similar to their SM values. Since the di- 
photon decay is dominated by the W^-boson contribution, the di-photon decay partial width is not 



10 




FIG. 2. M h>HjH ± vary with respect to M A for M~ t = 500 GeV, A t = -740 GeV, tan/3 = 11, fj, = 
2300 GeV. 

changed significantly from the SM r S M(-^ — > ll)- However, di-photon decay branching fraction 
BR(H — > 77) may still be enhanced due to decrease of H total width. At 125 GeV, H — >■ bb and 
H — > WW* H — > ZZ* dominate the H decay. Since H is mostly H u -\ike, H coupling to b is 
naturally suppressed. Given v u ~ v, gnzz and gnww are not significantly changed from the SM 
glf z and gf$ w . The partial widths of H ->• WW* and F ZZ* are indifferent from the SM 
values. With the reduction in H — > bb, the increase of H — > WW* and H — > ZZ* are inevitable. 
Therefore, light stau states in the spectrum can improve the di-photon behavior _R 77 and reduce 
the tension in increasing ZZ* or WW*. 

n 

Discussed in u£u, m the non-decoupling limit when H — >• bb still dominates the H decay, 
H — > t + t~ can be significantly enhanced. 

1 + A. 



*^H i+a^) J <14) 

where e = 1 + tana/ tan/3 with a < 0, A 6 is from the radiative correction in bottom Yukawa, 
r gg is the ratio in gluon fusion production of H which is order 1 in relatively large tan (3 and 
Mi > 500 GeV. With the radiative correction, H bb coupling is 



cos a 

gmb = n 

cos/3 



Af, / tana 



(15) 



1 + A h \ tan/3 

Similar to the story of [iA t in b — > s transition, A b also breaks Peccei-Quinn symmetry and R- 
symmetry at the same time. In this case, A b contains two i?-symmetry breaking pieces as gluino 
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mass Mg and A-term contribution. Our choice of fiA t < results in cancellation between the two 
contribution but the enhancement to R TT is still significant. Our results also confirm the finding 
in II lQM with many points of enhanced H — > tt decay. Figure [3] shows the correlation between 
BR(H —7- r + r~) and BR(H — > bb) in the survival points. On the other hand, we also find many 
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FIG. 3. BR(H -)• r+T~) in correlation with BR(H -)• bb). 



points with R TT < 1. One particularly interesting feature around non-decoupling limit is that 
H hh decay may open up and take significant portion of the H decay. In large parameter 
region, H — )■ hh decay partial width may completely dominate the decay of H once it opens up. 
Discussed in 114 ill , the tree level H ^ hh decay and one loop contributions may have different 
signs and severely cancel each otheiQ. There then exists a very fine tuned parameter region that 
the T(H — > hh) is at similar order as other decay and only takes about 50% of H decay. If 
H — )■ hh decay occurs, h can further decay into bb or rr, the search of H then fall into the 46,4r 
or 2b2r channels. The phenomenology of such channels have been widely studied in the context 
of NMSSM with h AA search 04QD. Studies of h ->■ AA in NMSSM shows that for M h ~ 
120 GeV, it requires the 14 TeV LHC with at least 100 fb _1 of data to claim discovery. Therefore, 
we argue the H — >■ hh decay is not constrained by any current direct search experimental data from 
LHC. In Fig 01 all the points R TT < 1 bare the same feature as BR(H — > hh) ~ 50%. Among 
these points, predictions on WW* and ZZ* are also slightly higher than the SM values but mostly 
within 1.5 which is consistent with the experimental data. The current search of H — > rr at 



The result is based on full one loop calculation in 14 ill and stability of the result may require higher order calculation. 
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ATLAS is still with large error bar and consistent with these large numbers of 2 a^. However, 
CMS collaboration has reported their latest data that exclude the SM tt rate by 1 a Ii42l1 . If one 
takes this seriously, most of our final survival parameter region will be cut away and only a few 
points that with significant H — > hh decay can survive. In addition, the H — >• bb are highly 
suppressed in these points and the predictions of these points agree with ATLAS central values of 
R in all channels very well. In principle, the choice of Ma can be extended to 0(200 GeV) in our 
study and the flavor bounds are less constrained for larger Ma- However, the larger Ma region 
corresponds to the enhanced R TT region. Only smaller Ma generates larger splitting between H 
and h which reduces R TT . Therefore, we only focus on the region Ma < 150 GeV. 

Besides the direct search via rr, LHC has put much stronger bounds on t — > bH + with H + — > 
t + u t comparing with Tevatron. The previous Tevatron upper bound of BR(t — > bH + ) is 5% while 
the latest ATLAS results become \%-5%. We plot the BR(t -> bH + ) with respect to M H ± by 
assuming BR(H + — > t + v t ) = 100% in FiglU It clearly shows that all the parameter points that 

0.020 
0.015 

x 

t 0.010 
0.005 
0.000 

120 130 140 150 160 

/W H+ (GeV) 

FIG. 4. BR(t -> bH + ) vs M H ± by assuming BR(# + — > t + v t ) = 100%. Red dots are parameter points 
that pass all our selection and constraints. 

pass our selections are below the search of light charged Higgs boson via top decay t bH + with 

H+ ->. T + V T . 

C. (7 X 7V 

Finally we discuss the last constraint for non-decoupling MSSM. Latest direct dark matter de- 
tection experiments XENON 100 have reached the level of sensitivity needed to detect neutralino 
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Data2011 7TeV 4.6fb _1 




dark matter over a substantial range of supersymmetric parameter space. These experiments at- 
tempt to detect weakly interacting (WIMP) dark matter particles through their elastic scattering 
with nuclei. Neutralinos can scatter with nuclei through both scalar (spin-independent) and axial- 
vector (spin-dependent) interactions. The experimental sensitivity to scalar couplings benefits 
from coherent scattering, which leads to cross sections and rates proportional to the square of 
the atomic mass of the target nuclei which is exactly being used for direct detection experiments. 
Consequently the spin-independent interactions are far more important than the spin-dependent in 
these experiments. In MSSM, the spin-independent interactions are mediated by the light Higgs 
bosons with cross section proportional to 

tan 2 /3 



(16) 



Figure \5\ has shown the spin-independent scattering between neutralino dark matter and the nuclei 
computed for XENON 100 setup by varying Ma for pure-bino of 200 GeV to illustrate the en- 
hancement feature due to light Higgs bosons and light top squarks. The calculation is done using 



micrOMEGAs 2.4 114311 . It clearly shows that the enhancement of such interaction in small Ma- In 



10 1 ' \ Mr = Mr = 500 GeV 



a 10 " 8 XENON 1 00(20 1 2) 



icr a 

Mr = Mr = 1000 GeV 



10 1 °130 140 150 160 170 180 

/W„(GeV) 



FIG. 5. Spin-independent scattering between neutralino dark matter and the nuclei computed for 
XENON100 by varying Ma- Red dashed line is in the case of light stop of 500 GeV and the blue dot- 
ted line corresponds to the stop mass of 1 TeV. 

addition, squark can induce neutralino-gluon scattering which can further enhance the scattering 



cross section 114411 . The two lines for different stop mass choices of 500 GeV and 1 TeV also indi- 



cates the enhancement of light stop in the neutralino-nuclei scattering. For most of our points with 
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500 GeV stop and Ma < 170 GeV, XENON100 bounds have put stringent constraints over the 
scenario. On the other hand, it is not clearly whether the current dark matter completely consists 
of supersymmetric neutralino. The bounds can also be easily evaded by adding new component of 
dark matter from non-supersymmetric origin. 



HI. CONCLUSIONS 



In this paper, we discuss the non-decoupling MSSM scenario where a light Higgs boson can 
evade the direct search experiments at LEP or Tevatron and the 125 GeV Higgs-like boson is 
explained as the heavy Higgs boson in the spectrum. The light Higgs boson may evade the di- 
rect search experiments at LEPII or Tevatron while the 125 GeV Higgs-like boson is identified as 
the heavy Higgs boson in the spectrum. Two direct consequences of the scenario are the flavor 
violation induced by the light charged scalar and the spin-independent scattering between neu- 
tralino and nuclei in dark matter direct detection experiments. With combined flavor constraints 
B — > Xg'j and B s — > ^ + ^" and direct constraints on Higgs properties, we find best fit scenarios 
with light stop of 0(500 GeV), negative A t around -750 GeV and large /i-term of 2-3 TeV How- 
ever, large parameter region in the survival space under all bounds may be further constrained by 
H — >■ rr if no excess of rr is confirmed at LHC. We only identify a small parameter region with 
significant H — >• hh decay that is consistent with all bounds and reduced rr decay. In addition, 
if current dark matter mostly consists of neutralino, direct detection experiments like XENON 100 
also puts stringent bound over this scenario with light Higgs bosons. The light stops which are 
required by flavor constraints can further enhance the scattering cross section. 



NOTE ADDED 



When completing our work, 1211.1955[hep-ph] [|45J] has appeared. The paper also studied 
similar region of non-decoupling MSSM and the results are in agreement with ours. We also 
include study on its enhancement of spin-independent neutralino-nuclei scattering. In addition, 
we find new parameter region which corresponds to reduce R TT due to H — > hh decay. 
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